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ABSTRACT 

It has been suggested that optically thin and geometrically thick accretion 
flows are responsible for the observed radio/X-ray luminosity relation of the X- 
ray bright galactic nuclei. If this is the case then central supermassive black 
hole masses can be estimated directly from measurements of the core radio lu- 
minosity and the X-ray luminosity, provided that properties of such accretion 
flows are known. Calculated ratios of the luminosities are presented in cases of 
the standard ADAF model and modified ADAF models, in which a truncation 
of inner parts of the flows and winds causing a reduction of the infalling matter 
are included. We compare the observed ratio of the luminosities with predic- 
tions from models of optically thin accretion flows. We also discuss the possible 
effects of the convection in ADAFs. We confirm that the supermassive black 
hole (SMBH) mass estimate is possible with the radio/X-ray luminosity relation 
due to ADAF models in the absence of a radio jet. We find that observational 
data are insufficient to distinguish the standard ADAF model from its modified 
models. However, the ADAF model with convection is inconsistent with obser- 
vations, unless microphysics parameters are to be substantially changed. High 
resolution radio observations are required to avoid the contamination of other 
components, such as, a jet component. Otherwise, the SMBH mass is inclined 
to be over-estimated. 

Subject headings: accretion disks - black hole physics - galaxies : nuclei - 
radiation mechanisms : nonthermal - radio continuum : galaxies - X-rays : 
galaxies 
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1. INTRODUCTION 

Since the early days of research on quasars and active galactic nuclei (AGNs), 
supermassive black holes (SMBHs) have been considered as the most likely power sources 
of the activity in these objects (Lynden-Bell 1969; Rees 1984). SMBHs at the centers of 
all galaxies are now recognised as ubiquitous, whose mass Msmbh is proportional to the 
spheroidal bulge mass of the host galaxy or the galactic bulge luminosity (Kormendy & 
Richstone 1995; Magorrian et al. 1998; Richstone et al. 1998). Evidence for the existence of 
SMBHs has been found in the center of our Galaxy (Eckart & Genzel 1997; Genzel et al. 
1997; Ghez et al. 1998) and in the weakly active galaxies NGC 1068, NGC 4258 (Miyoshi 
et al. 1995). Asymmetric Fe Ka emission in the X-ray spectra of AGNs (e.g., Tanaka et al. 
1995) may show a signature of SMBHs, but this remains somewhat speculative as Fe Ka 
reverberation signature have not yet been observed (e.g., Reynolds 2000). 

Searches for SMBHs are based on spatially resolved kinematics. In the case of AGNs, 
however, direct detection of nuclear SMBHs through stellar-dynamical methods has not 
been achieved due to technical difficulties. The reason is that bright AGNs are too bright 
to resolve the light from the surrounding stars and gas, and from the AGN itself on 
arcsecond and smaller angular scales. Reverberation mapping (Blandford & McKee 1982) 
avoids this problem. In this technique, the time delayed response of the emission lines to 
continuum variations is used to infer a size of the line-emitting region. We also have a 
velocity V, which is obtained from measurements of the Doppler widths of the variable line 
components. Combining these, a virial mass is estimated by Msmbh ~ V^R/G, where G is 
the gravitational constant. 

On the other hand, several authors have pointed out that reverberation mapping yields 
systemically smaller SMBH masses at a given bulge luminosity than do dynamical models 
of spatially resolved kinematics (e.g., Wandel, Perterson, & Malkan 1999). In the case of 
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normal galaxies, the SMBH mass appears to correlate with the galactic bulge luminosity, 
with the SMBH to spheriodal bulge mass ratio MsMBn/^buigc = 0.006 (Magorrian et al. 
1998; Richstone et al. 1998). In the case of AGNs, however, a significantly lower value of 
-^SMEnZ-^buigc has been found (Wandel 1999), indicating either a real difference between 
active and inactive galaxies, or that one or both of the methods of mass determination 
is somehow biased. Gebhardt et al. (2000) claim that masses derived from reverberation 
mapping are consistent with the recently discovered relation between the SMBH mass and 
the galaxy velocity dispersion derived from spatially resolved kinematics. 

Another possible way to infer the presence of SMBHs and to shed light on the physical 
condition is to examine spectral energy distribution over a wide range from the radio to 
the hard X-ray frequencies (e.g., Frank et al. 1992; Ho 1999). This emission spectrum is 
produced by an accreting matter, as the surrounding gas accretes onto the central SMBH. 
In such flows the core radio luminosity is low and dependent on the mass of the central 
SMBH (Ichimaru 1977; Rees et al. 1982; Narayan & Yi 1994, 1995a,b; Abramowicz et 
al. 1995). It has been suggested that advection-dominated accretion flows (ADAFs) are 
responsible for the observed radio and X-ray luminosities of some of the X-ray bright 
galactic nuclei (Fabian & Rees 1995; Di Matteo & Fabian 1997; Yi & Boughn 1998, 1999). 
In this case the central SMBH masses can be estimated directly from measurements of the 
X-ray and the core radio luminosities, provided that properties of such accretion flows are 
known. We demonstrate that it is relatively effective and consequently suitable to estimate 
SMBH masses of many candidates using a method we present in this paper. 

This paper begins with model descriptions for the standard ADAF model and modifled 
versions of the model in § 2. We deduce analytical expressions to describe the radio/X-ray 
luminosities for various models in § 2. We present the calculated ratios of the luminosities 
and compare them with the observational data in § 3. We discuss possible roles of the 
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convection in ADAFs and conclude in § 4. 

2. RADIO/X-RAY EMISSION FROM ACCRETING SMBHs 

When a mass accretion rate is below a critical rate, the radiative cooling is slower 
than the viscous heating. As a result of this, the dissipated accretion energy is inefficiently 
radiated away and advected inward to the central object with the accreted matter (see 
Narayan and Yi 1995b and references therein). The ions are heated by viscous dissipation 
at a rate of q^ per unit volume and the heating is parameterised by the viscosity parameter 
a. Since this flow is optically thin, the energy transfer from the ions to the electrons is 
inefficient so that the ions reach the virial temperature. The electrons cool via synchrotron, 
bremsstrahlung, and inverse Compton processes. Synchrotron radiation is responsible 
for the radio to submillimeter emission, bremsstrahlung emission and inverse- Compton 
scattering are submillimeter to X-ray emission. Detailed numerical calculations have been 
performed (Narayan et al. 1998), and the resulting spectra have been successfully applied 
to a number of extra galactic systems which are supposed to harbor SMBHs (e.g., Narayan 
et al. 1995; Lasota et al. 1996; Manmoto et al. 1997). 

In this study we apply analytical scaling laws for self-similar solutions (Narayan & Yi 
1994; Mahadevan 1997). We adopt the following dimensionless variables : mass of the SMBH 
m = M/Mq] radius from the SMBH r = R/Rg, where Rg = 2GM/c^ = 2.95 x lO^m cm; 
and mass accretion rate m = M/M^dd, where MEdd = LEdd/Vesc^ = 1-39 x lO^^m g s~^ (the 
Eddington accretion rate assuming r]eg = 0.1). We assume that the flows are spherically 
symmeric. The advection fraction / is determined such that the ion and electron energy 
balance equations are met. The quantities of interest are the volume-integrated quantities 
which are obtained by integrating the heating rate and each cooling rates throughout the 
volume of the flows. As canonical values in a model for ADAFs parameters are taken to be 
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rmin = 3, Vmax = 10^, a = 0.3, and P = 0.5 (see, e.g., Narayan & Yi 1995b). In this study, 
we assume that the electron temperature is constant for r < 10^, as suggested by Narayan 
& Yi (1995b). For given m,7h,a, (3, we may obtain the energy density spectrum with /, T^ 
determined as described below. 

2.1. Emission from 'Standard' ADAFs 

The ions are heated by the viscous dissipation and the electrons gain energy from the 
ions by Coulomb interactions alone. We neglect the possibility of heating the electrons, 
since we expect that the fraction of viscous energy transferred to the electrons is in the 
mass ratio of the electron to the ion, ~ 1/2000. 

The synchrotron photons are self-absorbed and give a blackbody radiation upto a 
critical frequency, Uc- For a given Te, the synchrotron spectrum L*^"'^ is given by 

Liy^^ = S3(siS2)'/'m'/'^4/'Tf/5z/2/5, (1) 
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where si = 1.42 x I0^a'^/^{1 - iSf'^c^ ^'^cj^, S2 = 1.19 x IQ~^^xm, and sg = 1.05 x 10 
xm = 'ii^ I^Vbdli ^b = eB/2'n-meC (Narayan & Yi 1995b; Mahadevan et al. 1996; Mahadevan 
1997). The radio luminosity Lf> at u is defined by uLiy. The highest radio frequency arises 
from the innermost radius of the accretion flows, Tmin ~ 3; z/p = SiS2m~^/'^m^/'^T^r^^l^ . At 
this peak frequency the peak radio luminosity is given by 

Ln = z/pLtf ^ = S3(sis,)3mV2^3/2y7^-7/4_ ^2) 

We obtain the total power due to synchrotron radiation from 

Psync ^ r L7^'du={hn)u,L:^;^. (3) 

Jo 
We set z/jnin = 0, since z/min is much smaller than z/p. 
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Bremsstrahlung emission is due to both electron-electron and electron-ion interactions. 
The total brenisstrahlung power is given by 

Pbrem = 4.74 X lO^^tt^^q ^ ln(r„„^/r^i„)F(^e)"^"^^ (4) 

and the spectrum due to bremsstrahlung emission is given by 

^brem ^ 2.29 X 10^' a-^c^Hn{r^a J rrmn)F {6, )mmX~' ^M'hiy / kT,) , (5) 

where 

9/3 

FiOe) = 4(^)1/2(1 + 1.781^,1-='^) + 1.73^f/2(i ^ 110^ ^ 02 _ i.25^f/2), 6, < 1, 

96' 

= (-^)[ln(1.123^e + 0.48) + 1.5] +2.3^e(lnl.l23^e + 1.28), ^e > 1- (6) 

27r 

In this study we neglect the Comptonization of bremsstrahlung emission, and consider 
the Comptonization of synchrotron emission alone. A contribution by Compton up-scattered 
synchrotron photons to the hard X-ray luminosity becomes important as m increases, while 
bremsstrahlung emission dominates the hard X-ray luminosity when m is substantially low. 

We approximate the Comptonized spectrum by assuming that all the synchrotron 
photons to be Comptonized have an initial frequency of Up. The maximum final frequency 
of the Comptonized photon is Uf = 3kTe/h, which corresponds to the average energy of the 
photon for saturated Comptonization in the Wien regime. On average, we assume that all 
the photons would see one half the total optical depth (Mahadevan 1997), which is written 
as Tes = Q.2a^^c^^mr~^/'^. The spectrum of the emerging photons at frequency v has the 
power-law shape 

L^''-^^L,,(-)-°^ (7) 

where ac = — Inres/lnA. The total Compton power is therefore given by 

3kTe/h 



J sync 



6.2 X 10^°Te 



Up 



1-Oc 

-1 



1 -ttc 

where A is the mean amphfication factor in a single scattering. 

For given m, m, a, /5, the total heating of the electrons should be balanced to the sum 
of the individual cooling, Q''^ = Psync + Phrem + Pcomp- The electron temperature is varied 
until this equality is satisfied. 

2.2. Emission from Truncated ADAFs and 'Windy' ADAFs 

ADAFs have the positive Bernoulli constant and, therefore, are susceptible to outflows 
(Narayan & Yi 1994). It might be the case that radio jets are present near the inner region 
of ADAFs. If so, the radio luminosity is dominated by jet emission. Instead, radio emission 
due to ADAFs at the high frequency will be supressed. In order to simulate this situation, 
we truncate ADAFs at a certain inner radius. We suppose that the truncation occurs at 
25rg. For instance, an observation of NGC 4486 (M87) indicates that the jet is formed in a 
smaller radius than r ^ SOr^ (Junor et al. 1999). This is slso about where is the maximum 
radius that 15 GHz radio emission could be generated for m = 10^, m = 10"^, Tg = 10^. 
Note that r = (1/60)-^/^ x (z//15 G}lz)'^/^{m/lO'^)'^/^{m/10-^y/^{TjlOy/^. We set r„i„ 
in equations in the previous subsection to 25rg instead of Sr^ in corresponding equations in 
this case. 

Since the gas in ADAF solutions is generically unbound, winds may carry away 
infalling matter (see, e.g., Blandford & Begelman 1999; Di Matteo et al. 1999; Quataert & 
Narayan 1999). We allow m to vary as m = rriout{T/ToutY (see, e.g., Blandford & Begelman 
1999). The mass loss causes a signiflcant effect on a model spectrum for a large value 
of p (Quataert & Narayan 1999). Bremsstrahlung emission decreases with increasing p, 
synchrotron emission decreases strongly with increasing p due both to the lower density 
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and to the lower Tg near r ~ 3. Compton emission decreases with increasing p even more 
strongly than the other two emissions since it depends both on UpL,^ and ac- A wind model 
is incompatible with observation, if microphysics parameters we adopt in the standard 
ADAF model are more or less correct. This conclusion, however, is inconclusive since there 
are qualitative degeneracies between the wind parameter, p, and microphysics parameters, 
a,f3, particularly the viscous heating paramter of the electrons, 6 (Quataert & Narayan 
1999). As rh is replaced by rhout{f^/routy equations for the 'Standard ADAFs' are to be 
modified as described below. The total heating for ions Q~^ is rewritten as 

Q+ = 9.39xl03«iy^C3mm,„,r;j;-l^(C-i-r:T^), (9) 

where p < 1. Similarly, Q^'^ is given by 

Q- = 1.2 X lO''9{e.)a-\'csPmml,r^r^^{r'j:a.' ' r'Zn) (10) 

where p ^ 1/2. We assume that xm and Te are constants with r and m for simplicity 
as a first order approximation. Given xm, the cutoff frequency at each radius is given as 
Uc = SiS2m~^^'^mJu^fr~^/ j'2^-(2p-5)/4_ rjj_^g synchrotron spectrum is given by 

Lr^ = S3(siS2)^mi^<#r,#Te^z/t^. (11) 

The peak frequency and the radio luminosity at the peak frequency are correspondingly 
reduced so that we have 

lypLtl^^ = Ss{s,S2)'m'/'m'J,yJ;TjrS . (12) 

The total synchrotron power is also rewritten by 

Psync ^ j^ LT^du={^—^)UpLiy;^. (13) 

For bremsstrahlung emission the total power and the spectrum are respectively given by 

1 

% 

^,rem ^ 2.29 X lO^V" V^mLr;„ti^(^e)T-i exp(-/^z//A;Te) ^(i - r^^J]. (15) 



\rem = 4.74 X 10^'' a'' cf mml,r-;^r F (9,) [-{r'j:^^ - r'^J], (14) 
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Note that as the exponent p approaches to 0, [^(^max~''"mm)] becomes ^^{rmax/^min)- For the 



Comptonization of the synchrotron photons we take t^s as t^s = Q-2a^^ Ci^moutrouti''"rnax — 
'''rrdn )• "^^^ Compton spectrum and the total Compton power are respectively given by 

L^"^P ~ ^.,(-)-"% (16) 



-' comp , 

1 — a 



i^Mr' rr6.2xioioT,^'-"^ 



Up 



-1 . (17) 



3. BLACK HOLE MASSES AND RADIO/X-RAY LUMINOSITIES 

Figure 1 shows a plot of the ratio of the 15 GHz radio luminosity to the 2-10 keV 
X-ray luminosity versus the X-ray luminosity. The solid lines are standard ADAF model 
predictions for SMBH masses of 10^ — IO^Mq, the dotted lines are ADAF model predictions 
with truncations at 25rg. The short dashed lines and the long dashed lines are ADAF model 
predicitons with winds, where p = 0.4 and p = 0.99, respectively. When p = 1, results may 
represent a case where the convection is present in ADAFs (Quataert & Narayan 1999). 
For each curve, m or rhout for windy models varies from 10~^ (upper left corner) to 10~^'^ 
(lower right corner). The open circles denote the observed core radio luminosity at 15 
GHz, and the filled circles the spatially resolved core radio luminosity which are converted 
from 5 GHz to 15 GHz using the z/^/^ power law. When the accretion rate m is large the 
change of the properties of the flows is considerable in the ratios. As long as the canonical 
model parameters remain unchanged, the ADAF models with large p, which represent the 
convection, seem to be ruled out by the observational data. In Figure 2, the core radio 
luminosity is shown as a function of the SMBH mass, with SMBH candidates whose mass 
estimates are available. From top to bottom, for each set of lines, the mass accretion rate 
corresponds to 10^^ and 10~^. The open circles indicate the observed core radio luminosity 
at 15 GHz, and the filled circles the converted core radio luminosity from 5 GHz to 15 GHz, 
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as in Figure 1. Although the radio flux data is in high angular resolution, there appears a 
somewhat wide range of radio luminosities, which is likely the result of radio jets or other 
components of various frequencies rather than pure ADAF emission. In Figure 3, the X-ray 
luminosity for the 2-10 keV band is shown as a funciton of the SMBH mass. 

Radio flux values at 15 GHz of Sgr A*, NGC 1068 are quoted from Kormendy & 
Richstone (1995); Gallimore et al. (1996), NGC 3079, NGC 3628, NGC 4151, and NGC 
4388 are quoted from Carral et al. (1990), NGC 3031 (M81), NGC 4258, NGC 4736 (M94), 
and NGC 5194 (M51) are quoted from Turner & Ho (1994), respectively. Radio flux data 
at 5 GHz of NGC 1365 and NGC 3310 are quoted from Saikia et al. (1994), and NGC 224 
(M31), NGC 3377, NGC 4374 (M84), NGC 4486 (M87), and NGC 4594 (M104) are quoted 
from Franceschini et al. (1998), respectively. 

X-ray fluxes (2-10 keV) of Sgr A*, NGC 224, NGC 1068, NGC 1365, NGC 3031, 
NGC 3079, NGC 3310, NGC 3628, NGC 4151, NGC 4258, NGC 4374, NGC 4388, NGC 
4594, NGC 4736, and NGC 5194 are quoted from Koyama et al. (1996); Trinchieri et al. 
(1999); Turner et al. (1997); lyomoto et al. (1997); Ishisaki et al. (1996); Cagnoni et al. 
(1998); Ptak et al. (1999); Ptak et al. (1998); Nandra et al. (1997); Makishima et al. (1994); 
Colbert & Mushotzky (1999); Iwasawa et al. (1997), Forster et al. (1999); Ptak et al. 
(1998); Roberts et al. (1999); Terashima et al. (1998), respectively. Quoted X-ray fluxes 
of NGC 4486 and NGC 3377 should be considered as upper limits (Reynolds et al. 1999; 
Pellegrini 1999). Note that the observed X-ray flux of NGC 3377 in the energy band of 0.2 
- 4 keV is converted to the value in the energy band of 2 - 10 keV, assuming a power-law 
with the canonical spectral index for Seyfert 1 galaxies, i.e., F = 1.7. High resolution X-ray 
images have revealed discrete X-ray sources for NGC 224 (Trinchieri et al. 1999), NGC 
3031 (Ishisaki et al. 1996), NGC 4736 (Roberts et al. 1999). For these sources we adopt 
core or bulge X-ray fluxes. 
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The mass estimates of Sgr A*, NGC 1068, NGC 3031 (M81), NGC 3079, NGC 4374 
(M84), 4594 (M104) come from Eckart & Genzel (1997); Greenhill et al. (1996); Ho (1999); 
Ferrarese & Merritt (2000); Ho (1999); Ferrarese & Merritt (2000), respectively. NGC 224 
(M31) is one of strongest SMBH candidates because of the fast rotation and the large 
velocity dispersion (Richstone et al. 1990). The estimates of SMBH mass come from 
Dressier & Richstone (1988), Bacon et al. (1994). NGC 224 seems to harbour a double 
nucleus (Lauer et al. 1993; Bacon et al. 1994). The mass estimate of NGC 3377 comes from 
Kormendy & Richstone (1995) and Magorrian et al. (1998). Besides M32, NGC 3377 is 
the second elliptical galaxy with stellar dynamical evidence for SMBHs. The mass of NGC 
4151 have been measured by reverberation mapping method (Wandel et al. 1999; Gebhardt 
et al. 2000). The mass estimate of NGC 4258 comes from Miyoshi et al. (1995). NGC 
4258 belongs to LINERs (e.g., Osterbrock 1989), which is one of successful applications of 
ADAFs (Lasota et al. 1996). The mass estimates of NGC 4486 (M87) come from Ford et 
al. (1994), Reynolds et al. (1996), Magorrian et al. (1998), and Ho (1999). NGC 4486 is 
a radio-loud elliptical galaxy with an optical jet which is approximately perpendicular to a 
disk of ionized gas. 

4. DISCUSSION AND CONCLUSION 

One of important predictions of the ADAF model is the radio/X-ray luminosity relation 
which can be directly used to estimate central SMBH masses. We demonstrate with an 
improved statistical significance that if ADAFs are confirmed by observations with a high 
angular resolution, the inherent radio/X-ray luminosity relation provides a direct estimate 
of the central SMBH mass. Several nearby extra galaxies have consistent ratios of radio 
to X-ray luminosities with the ADAF predictions for an estimated mass of the central 
SMBH. We plot ADAF model predicitons for the radio/X-ray luminosities in terms of the 
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SMBH mass in Figure 1. The observations for several extra galactic objects are consistent 
with the predictions of ADAF models, as seen in Figures 2 and 3. Although the mass 
estimates of NGC 4374, NGC 4486, NGC 4594 appear more consistent with the predictions 
with truncated ADAF models or windy ADAF models, there is no obvious signs of the 
truncation or the wind in ADAFs. It is, however, interesting to note that the models with 
the convection are exclusive unless the microphysics in the standard ADAF model has to 
be modified significantly. 

The predicted mass of NGC 224 (M31) is smaller than the estimated mass at least 
by two orders of magnitude. NGC 224 has an extrememly low core radio luminosity, 
which is even lower than predicted with the mass accretion rate that the X-ray luminosity 
implies on the basis of an ADAF model. NGC 224 seems to harbor a double necleus 
(Lauer et al. 1993; Bacon et al. 1994). It is unclear whether ADAF is viable under such 
circumstance. Predicted masses of NGC 1068, NGC 3031 (M81), NGC 3079, NGC 4151 
are ~ 10® to ~ a few 10^ Mq. On the other hand, reported mass estimates of SMBHs 
of these galaxies are less than ~ 10'^. Masses of these sources are over-estimated due to 
high radio luminosities resulting from strong jet-related activities. Moreover, uncertainties 
in X-ray flux measurements and intrinsic variations could cause errors in our predictions. 
NGC 1068 is famous in that it has a strong compact radio source, which results in a much 
higher radio/X-ray luminosity ratio than other Seyferts (Gallimore et al. 1996). The radio 
core of NGC 3031 (M81) is highly variable on many timescale (Ho 1999). Moreover, NGC 
3031 is substantially brighter in X-rays relative to the UV than in luminous AGNs (Ho 
1999), and its core luminosity in the 2-10 keV band varies by a factor of ~ 2 in a timescale 
of years (Ishisaki et al. 1996). The unresolved radio map of NGC 3079 at 15 GHz shows 
an intriguing feature, which appears just like a jet (Carral et al. 1990). The high radio 
luminosity may be due to an unresolved small-scale jet. The radio morphology of NGC 
4151 suggests that the emitting material has been energized from the nucleus and that the 
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jet interacts with the ambient medium (Carral et al. 1990). 
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Fig. 1. — A plot of the ratio of the radio luminosity to 2-10 keV X-ray luminosity versus the 
X-ray luminosity. The sources are spatially resolved, and the core luminosities are adopted in 
order to avoid possible contaminations due to other components such as jets. SMBH masses 
are denoted by the log scale at the top of each curve. See text for detailed discussions and 
references. 

Fig. 2. — The core radio luminosity is shown as a function of the SMBH mass, with SMBHs 
whose masses are available. Mass accretion rates are denoted by the log scale at each line. 
Lines correspond to same models as in Figure 1 and symbols represent same observational 
data points as in Figure 1. 

Fig. 3. — Similar plot as Figure 2. But the X-ray luminosity for the 2-10 keV band is shown, 
with the mass-known SMBHs. Lines correspond to same models as in Figure 1. 
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